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ABSTRACT 

An  adaptive  thresholding  method  is  presented  for  optimum  detection  for  optical  receivers  with  large  multiplicative 
noise.  The  technique  uses  low  frequency  sampling  of  the  detected  current  that  enables  calculation  of  the  bit  means  and 
variances  and  estimation  of  the  optimum  detection  threshold.  The  regime  in  which  this  holds  is  when  the  sampling 
frequency  is  lower  than  the  bit  rate  but  higher  than  atmospheric  turbulence  frequency  content.  Simulations  are  done 
with  data  obtained  from  the  NRL  Chesapeake  Bay  Lasercomm  Testbed.  The  results  of  simulations  comparing  BER 
performance  versus  sample  rate  and  parameter  estimation  error  will  be  presented.  If  the  system  parameters  are 
characterized  in  advance  with  reasonable  accuracy,  the  BER  obtained  will  typically  be  an  order  of  magnitude 
improvement  over  the  equal  variance  threshold  (depending  on  the  signal  to  noise  ratio). 
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1.  INTRODUCTION 

It  is  well  known  that  the  optimum  detection  threshold  for  optical  receivers  with  large  multiplicative  noise  components 
(i.e.  signal  dependant  noise)  is  not  simply  the  average  value  of  the  high  and  low  bit  currents  (or  zero  for  ac  coupled 
detectors).  The  optimum  detection  threshold  in  this  type  of  system  is  derived  from  a  Bayes'  likelihood  ratio  test  (LRT) 
and  is  a  function  of  the  bit  level  means  and  variances;  therefore,  it  will  not  be  constant  under  varying  transmission 
conditions  [1,2].  Optical  detectors  exhibiting  this  characteristic  are  avalanche  photodiodes  (APD’s)  and 
photomultipliers  as  well  as  PIN  photodiodes  coupled  with  optical  preamplifiers  [2, 3, 4, 5].  The  derivation  of  the 
detection  threshold  requires  a  priori  knowledge  of  the  mean  and  variance  of  the  bit  levels.  As  stated  above,  in  free-space 
optical  communication  links  this  will  be  problematic.  An  adaptive  method  that  can  track  the  changes  in  the  mean  and 
variance  of  the  signal  bits,  and  update  the  detection  threshold  is  needed  in  order  to  implement  the  multiplicative  noise 
Bayesian  LRT  detection  threshold  [3,6]. 

In  free-space  optical  communication  (FSO)  systems,  it  is  usually  desirable  to  have  the  sensitivity  of  the  detector  as  high 
as  possible  to  reduce  the  required  laser  power  for  a 

specified  link  margin.  However,  the  optical  receivers  used  for  high  speed  optical  communication  systems  that  operate 
closest  to  the  quantum  limit  are  the  PIN  diode  coupled  to  an  optical  fiber  preamplifier  and  an  APD.  Both  of  these 
devices  have  large  multiplicative  noise  components  that  cannot  be  ignored. 
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In  previous  work,  several  adaptive  predictor  methods  have  been  developed  which  estimate  both  the  mean  and  the 
variance  of  the  bit  levels  well  enough  to  produce  bit  error  rate  performance  at  the  theoretical  limit  [7,8].  These  methods 
include  Kalman  filtering  for  single  or  multiple  samples  per  bit,  least  mean  squares  (LMS)  adaptive  algorithm,  and  a 
modified  sequential  regression  (SER)  adaptive  algorithm.  These  adaptive  signal  processing  algorithms  require  high 
speed  sampling  equipment  and  high  speed  computers  for  implementation.  Since  this  type  of  equipment  can  be  costly 
and  difficult  to  implement,  it  is  desired  to  find  a  method  of  implementing  adaptive  thresholding  that  is  simple  and 
reasonably  cost  effective.  The  purpose  of  the  work  discussed  in  this  paper  is  to  develop  an  inexpensive  implementation 
method  for  a  lower  bit  rate  optical  link  such  as  an  asymmetric  100  Mbps  link  based  on  a  multiple  quantum  well  (MQW) 
modulating  retroreflector  device  [9].  Since  these  lightweight  devices  are  possible  candidates  for  portable  links,  an  APD 
would  increase  the  sensitivity  of  the  interrogator  without  adding  weight,  while  reducing  the  required  transmitter  laser 
power. 


2.  BACKGROUND 


As  stated  above,  the  optimum  detection  threshold  for  an  optical  detector  with  signal  dependant  noise  that  cannot  be 
ignored  can  be  derived  from  a  Bayes’  Likelihood  Ratio  Test  [7,8].  The  expression  for  the  optimum  threshold  is 
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in  which  q0  and  |li  are  the  mean  signal  currents  of  the  low  and  high  bits  respectively  and  the  G02  and  Gi2  are  the 
respective  variances.  If  the  signal  dependant  noise  is  negligible  as  it  is  in  a  standard  PIN  diode  (with  no  optical 
amplifier),  g02  and  Gi2  are  approximately  equal  and  the  Eqn.  (1)  will  reduce  to  the  average  value  of  q0  and  This 
approximation  is  the  equal  variance  threshold  (EVT). 


The  APD  has  an  increased  shot  noise  due  to  the  avalanche  gain  process.  The  APD  shot  noise  current  variance  is 
modeled  as: 

=  1  eC;-F(Gm\RPrB,  (2) 

where  Gm  is  the  avalanche  multiplicative  gain,  Pr  is  the  received  optical  power,  Be  is  the  electronic  bandwidth,  R  is  the 

detector  responsivity  at  a  gain  of  1 ,  and  F( Gm)  is  the  excess  noise  factor  which  is  a  function  of  the  avalanche  gain.  The 
APD  shot  noise  is  typically  modeled  as  a  zero  mean,  Gaussian  process  at  the  signal  levels  encountered  in  optical 
communication  systems.  The  excess  noise  factor  is  given  by: 

F  (GJ  =  •  G,  +  ( 1  -  kA )  ■  (2  -  yG J  (3) 

where  kA  is  the  ionization  coefficient  which  is  a  property  of  the  semiconductor  material  used  in  the  APD  and  has  a  range 
from  0  to  1.  In  silicon  APD's,  which  can  be  used  up  to  an  optical  wavelength  of  ~  llOOnm,  the  ionization  coefficient 
can  be  kept  fairly  small  (typically  «  1).  However,  for  optical  wavelengths  above  -  1100  nm,  such  as  1330  nm  and 
1550  nm  which  are  of  particular  interest  for  freespace  links,  other  semiconductor  materials,  such  as  InGaAs,  are  used. 
InGaAs,  APD's  have  an  ionization  coefficient  much  higher  than  that  of  silicon,  typically  in  the  vicinity  of  ~  0.7  [10,11]. 

Therefore,  since  the  signal  dependent  shot  noise  in  APD's  is  multiplied  by  an  additional  factor  of  Gm 2  •  F  (Gw)  ,  the 
shot  noise  in  these  devices  will  dominate  over  thermal  noise  and  cannot  be  neglected. 


In  the  PIN  diode-optical  preamp  receiver,  there  will  be  additional  terms  in  the  current  noise  variance  due  to  the  beating 
of  the  spontaneous  noise  power  produced  in  the  optical  amplifier  with  itself  and  with  the  received  optical  signal.  These 
noise  terms  can  also  be  adequately  modeled  as  Gaussian  processes  and,  at  typical  optical  amplifier  gains,  will  be  the 
dominant  noise  terms  over  both  thermal  noise  and  shot  noise.  The  signal- spontaneous  beat  noise  produces  a  noise 
current  variance  of: 

<7  ,2-4  R2  G  P  P  (G  -l)  B  (4) 

sig-spont  o  r  n  \  o  )  e  v  ' 
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where  G0  is  the  optical  amplifier  gain  and  Pn  =  n  •  h  •  Vc  ,  in  which  /z  is  Planck's  constant,  V c  is  the  optical  frequency, 

and  nsp  is  the  spontaneous  emission  factor  (typically  in  the  range  of  2  to  5)  [10].  The  spontaneous-spontaneous  beat 
noise  produces  a  noise  current  variance  of: 

=  2  R2  ■  [p„  •  (G„  - 1)]!  •  (2 .  -  B, )  ■  B,  (5) 

in  which  B0  is  the  optical  bandwidth  of  the  system,  usually  determined  by  optical  bandpass  filters.  In  practice,  to  reduce 
the  system  noise  and  improve  performance,  the  spontaneous- spontaneous  beat  noise  term  can  usually  be  reduced  to  a 
very  small  level  with  optical  filters  so  that  the  dominant  noise  term  is  always  the  signal-spontaneous  beat  noise.  As  can 
be  seen  above,  this  dominant  noise  term  is  also  proportional  to  the  received  optical  power. 

As  shown  above,  the  most  sensitive  detector  configurations  for  freespace  optical  communication  systems  have  large 
multiplicative  (signal  dependent)  noise  variances  which  cannot  be  ignored.  Although  the  PIN-diode-optical-preamp 
receiver  has  a  sensitivity  closest  to  the  quantum  limit,  the  APD  receiver  has  a  much  wider  use  in  FSO.  Currently, 
optical  preamplifiers  are  constructed  with  single-mode  optical  fiber.  In  freespace  optical  systems,  turbulence-induced 
spot  motion  makes  coupling  of  received  light  into  a  single-mode  fiber  extremely  difficult  and  the  sensitivity  gain  is 
canceled  by  the  coupling.  Therefore,  the  analysis  and  simulations  done  for  this  paper  concentrates  on  systems  using 
APD  detectors. 

Atmospheric  turbulence  can  cause  the  average  received  power  to  fluctuate  with  power  spectrum  components  in  the 
kilohertz  regime.  Therefore,  the  received  bit  signal  currents  and  their  variances  fluctuate  as  well.  In  order  to  maintain 
the  receiver  detection  threshold  at  near-optimum,  these  bit-level  means  and  variances  must  be  tracked  and  estimated.  A 
block  diagram  of  a  proposed  solution  is  shown  below  in  Figure  1 . 


Previous  work  has  involved  theoretical  comparisons  of  bit  error  rates  obtainable  with  the  optimum  Bayesian  LRT 
detection  threshold  and  the  equal  variance  threshold.  An  example  of  this  sort  of  analysis  is  shown  in  Figure  2  in  which 
it  can  be  seen  that  more  than  an  order  of  magnitude  improvement  can  be  obtained  with  the  adaptive  LRT  threshold  [7,8]. 
It  has  also  been  shown  that  use  of  an  LRT  threshold,  established  at  some  fixed  value,  and  not  updated  over  time,  will 
actually  perform  worse  than  an  EVT  threshold  [7,8]. 


Proc.  ofSPIE  Vol.  5160  357 


These  past  efforts  [7,8]  have  concentrated  on  development  of  adaptive  predictor  algorithms  which  can  predict  the  mean 
and  variance  values  and  maintain  the  detection  threshold  near  optimum.  Adaptive  algorithms  have  been  developed 
based  on  Kalman  Filters,  Least  Mean  Squares  adaptive  predictors,  and  a  Modified  Sequential  Regression  adaptive 
algorithm.  It  has  been  demonstrated  with  simulation  studies  that  these  adaptive  predictors  can  yield  BER  performance 
at  the  theoretical  limit,  with  order-of-magnitude  improvement  over  the  EVT  threshold  performance.  Figure  3  shows  a 
comparison  of  the  three  adaptive  predictors’  BER  performance  to  the  theoretical  limit.  Figure  4  is  a  comparison  of  the 
Kalman  Filter  adaptive  predictor  BER  performance  to  that  of  the  EVT  threshold.  The  results  were  generated  with 
simulated  bit  data  superimposed  on  actual  power  fade  rate  data  taken  at  the  NRL  Free-space  Lasercomm  Test  Facility  at 
Chesapeake  Beach,  MD  [7,8]. 


358  Proc.  of  SPIE  Vol.  5160 


Adaptive  Predictor  BER  Comparison 


•  Theoretical 
^^“0  Errors 

Kalman  Filter 


LMS 

^I^SER 


Figure  3:  Comparison  of  Theory,  Linear  State  Model  Kalman  filter,  LMS  predictor  and  Modified  Sequential 
Regression  predictor  bit  error  rate  performance;  averaged  over  2msec  time  intervals. 


Kalman  Filter  /  EVT  :  BER  Comparison 


Figure  4:  Comparison  of  Kalman  filter  BER  performance  with  the  EVT. 
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3.  LOW  FREQUENCY  SAMPLING 


Implementation  of  the  adaptive  predictors  described  above  requires  high  speed  sampling  and  high  speed  computing 
ability.  These  processing  requirements  can  be  very  costly  and  difficult  to  realize. 

It  is  desired  to  find  a  way  to  implement  adaptive  thresholding  that  is  cheaper  and  easier.  If  the  free-space  optical  link 
components  can  be  well  characterized  in  advance,  it  is  possible  to  estimate  the  optimum  detection  threshold  by 
estimating  the  average  received  detector  current  before  AC  coupling.  For  a  communication  link  using  an  APD,  the 
parameters  which  must  be  known  or  estimated  in  advance  are:  transmitted  contrast  ratio,  non-multiplicative  circuit  noise 
variance,  APD  gain,  APD  ionization  ratio,  APD  excess  noise  factor,  and  APD  bandwidth.  These  parameters  are 
accessible  from  manufacturer’s  specifications  or  measurement.  With  these  parameters  known  or  estimated,  and  with  a 
continuous  measurement  of  the  receiver  detector  current,  the  high  and  low  bit  mean  and  variance  levels  can  be 
calculated.  Consequently,  calculation  (estimation)  of  the  optimum  detection  threshold  is  enabled. 

A  testbed  is  being  developed  at  the  Naval  Research  Laboratory  to  test  the  feasibility  of  this  method  in  the  laboratory. 
The  average  receiver  current  will  be  estimated  by  low  frequency  sampling  of  the  received  signal  before  AC  coupling 
with  a  commercial,  off-the-shelf  (COTS)  Analog-to-Digital  (A/D)  card  from  National  Instruments.  A  block  diagram  is 
shown  in  Figure  5.  System  parameters  will  be  characterized  in  advance  and  the  estimated  optimum  detection  threshold 
will  be  updated  with  a  COTS  D/A  card  and  a  COTS  Data  and  Clock  Recovery  chip.  Data  and  Clock  Recovery  chips 
may  be  obtained  that  have  a  control  voltage  input  for  adjusting  the  detection  threshold  used  in  the  data  recovery. 
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Figure  5:  Block  diagram  of  low  frequency  sampling  testbed  being  developed  at  NRL. 


4.  SIMULATION  RESULTS 


As  previously  stated,  studies  were  performed  with  data  consisting  of  a  random  bit  pattern  imposed  on  power  fade  rate 
data  taken  at  the  NRL  Lasercomm  Test  Facility.  The  average  receiver  current  was  estimated  at  sampling  rates  of  1,  2,  5, 
10,  50,  100,  and  500  kHz  using  the  simulated  input  data.  The  estimated  receiver  current  was  then  used  to  calculate  the 
estimate  of  the  optimum  detection  threshold.  The  estimated  thresholds  were  used  to  detect  the  bits  on  the  AC  coupled 
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data.  The  detected  bits  were  compared  to  the  true  transmitted  bit  sequence  and  the  BER  was  compared  to  that  obtained 
with  the  true  optimum  detection  threshold.  These  studies  of  BER  as  a  function  of  sampling  frequency  were  performed 
for  two  average  signal-to-noise  ratios  (SNR)  of  17  and  43.  The  different  SNR  values  were  obtained  by  scaling  the 
original  digitized  fade  rate  data  before  generating  the  simulation  data.  The  simulation  data  (4  million  bits)  was 
generated  by  the  process  depicted  in  Figure  6.  The  parameters  used  in  the  data  generation  process  correspond  to  COTS 
components;  specifically:  an  OCP  STX-12  transmitter  module  operating  at  1550nm  with  a  contrast  ratio  of  lOdB  and  a 
Lucent  Technologies  126C  InGaAs  APD. 


Figure  7  contains  a  graph  of  the  true  optimum  detection  threshold  for  the  case  of  average  SNR  of  17.  In  the  graph,  it  is 
shown  that  the  detector  threshold  value  is  in  the  microamps  regime.  Although  the  simulations  were  done  using  detector 
current,  in  reality,  the  detector  current  will  be  multiplied  by  the  trans -impedance  gain  of  the  trans -impedance  amplifier 
(TIA).  This  factor  has  no  effect  on  the  simulation  outcome  (other  than  a  slight  decrease  in  SNR)  since  the  TIA  gain  is 
simply  a  multiplier  of  the  current  and  threshold  for  implementation.  Figure  7  is  the  true  optimum  detection  threshold  for 
the  data  set  being  used  for  the  average  SNR  =17  case  and  should  be  compared  with  the  estimates  obtained  by  low 
frequency  sampling  in  figures  9  and  10. 

Figure  8  is  an  example  of  the  AC  coupled  data  that  is  used  in  the  simulation  experiment  at  an  average  SNR  of  17.  The 
plot  in  Figure  8  plots  every  tenth  data  sample  in  the  4  million  generated.  The  power  fluctuations  caused  by  turbulence 
are  apparent  as  are  the  effects  of  the  multiplicative  noise;  i.e.  the  “one”  bits  are  more  noisy  than  the  “zero”  bits  and  the 
plot  is  not  symmetric  about  the  x-axis. 
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_ Figure  8:  AC  coupled  received  data;  SNR  ave  =17. _ 


Figures  9  and  10  are  examples  of  the  estimated  receiver  current  and  estimated  optimum  detection  threshold  for  SNR=17, 
for  two  different  sample  rates,  5  and  500  kHz.  Estimated  detection  thresholds  were  obtained  for  a  range  of  sampling 
frequencies  from  1  kHz  to  500  kHz.  These  were  determined  from  the  estimates  of  the  receiver  current  and  the 
predetermined  or  estimated  values  of  the  transmitter  and  receiver  parameters  such  as  contrast  ratio  and  APD  gain. 
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Figure  9:  Estimated  receiver  current  and  estimated  detection  threshold;  Sample  frequency  =  5  kHz;  SNR 

ave  =  17. 
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Figure  10:  Estimated  receiver  current  and  estimated  detection  threshold;  Sample  frequency  =  500  kHz; 

SNR  ave  =17. 
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To  obtain  the  data  shown  in  Figures  11  and  12,  the  estimated  detection  thresholds  for  the  various  sampling  frequencies 
were  used  to  detect  the  bits  in  the  AC  coupled  data.  These  detected  bits  were  then  compared  to  the  true  transmitted  bit 
stream  to  determine  average  bit  error  rate.  The  data  was  scaled  to  obtain  SNR  values  that  would  result  in  average  bit 
error  rates  greater  than  1  x  10‘4.  This  is  the  region  in  figure  2  in  which  the  theoretical  curves  for  the  optimum  threshold 
and  for  the  EVT  threshold  are  fairly  close  together  so  that  small  changes  in  BER  performance  are  obvious.  The  lower 
straight  line  in  these  figures  is  the  BER  obtained  by  using  the  true  optimum  detection  threshold.  The  BER  obtained  by 
using  the  EVT  threshold  is  listed  in  the  upper  right  corner  of  the  graphs.  As  can  be  seen  in  the  figures,  if  the  sampling 
frequency  of  the  received  current  monitor  is  10  kHz  or  higher,  the  BER  performance  will  be  very  close  to  that  obtained 
with  the  true  optimum  detection  threshold  and  depending  on  the  SNR,  an  order  of  magnitude  improvement  would  be 
obtained  compared  to  the  EVT  threshold. 
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Since  certain  system  parameters  discussed  above  are  required  to  be  known  a  priori  in  order  to  calculate  the  detection 
threshold,  it  is  desirable  to  have  an  understanding  of  the  accuracy  with  which  these  parameters  need  to  be  known.  Since 
APD  gain  and  the  contrast  ratio  of  the  transmitter  have  a  large  effect  both  on  the  mean  and  the  variance  of  the  bit  levels, 
simulation  studies  were  performed  in  which  these  parameters  were  intentionally  changed  from  the  true  value  to 
determine  the  effect  on  the  detection  threshold  and  BER.  The  studies  were  done  for  two  different  average  SNR  values 
(12  and  59)  and  for  parameter  estimation  errors  of  +/-50%.  The  sampling  frequency  used  for  these  simulations  was  50 
kHz  and  the  total  number  of  bits  used  was  500,000. 

Figure  13  shows  the  results  for  these  parameters.  The  side-by-side  graphs  contain  the  same  data  but  with  different 
scales  so  that  a  comparison  to  the  EVT-threshold  BER  could  be  done  in  the  graph  on  the  right.  The  higher  average  SNR 
case  (SNR  =  59)  has  a  two  level  graph  that  is  the  difference  between  one  error  and  two  errors  out  of  500,000  bits.  This 
indicates  that  in  the  high  SNR  case,  the  BER  is  not  very  sensitive  to  errors  in  the  detection  threshold  due  to  gain 
estimation  errors.  Figure  14  is  the  simulation  results  for  errors  in  the  estimate  of  the  contrast  ratio.  As  in  Figure  13,  the 
upper  two  graphs  are  for  the  lower  average  SNR  of  12  and  the  lower  two  graphs  are  for  the  average  SNR  of  59.  It  is 
shown  that  an  overestimate  of  the  contrast  ratio  can  result  in  a  large  increase  in  the  BER,  as  one  would  expect. 

It  should  be  noted  that  in  the  analysis,  the  received  current  estimates  were  obtained  by  averaging  over  the  integration 
time  of  the  A/D  card.  It  was  assumed  that  the  integration  time  was  equal  to  the  inverse  of  the  sample  frequency.  An 
attempt  was  made  to  use  an  LMS  linear  transversal  filter  as  a  predictor  with  the  averaged  current  values  as  inputs  to 
improve  the  estimation  of  the  detection  threshold.  However,  this  resulted  in  little  or  no  improvement  in  the  BER,  so  the 
measured  average  current  value  was  used  as  the  prediction  of  received  current. 
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BER  vs.  APD  Gain  Estimate  Error,  SNR  ave  =  12 
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Figure  13:  Effects  of  APD  gain  estimate  error. 


366  Proc.  of  SPIE  Vol.  5160 


BER  vs.  Contrast  Ratio  Estimate  Error 
SNR  ave  =  12 
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Figure  14:  Effects  of  contrast  ratio  estimate  error. 


5.  CONCLUSIONS 

With  low  frequency  sampling  of  the  received  detector  current,  an  estimation  of  the  detection  threshold  can  be  obtained 
which  is  usefully  comparable  to  that  obtained  with  the  true  optimum  detection  threshold.  If  the  system  parameters  are 
characterized  in  advance  with  reasonable  accuracy  (+  /-  10%),  the  BER  obtained  will  typically  be  an  order  of  magnitude 
improvement  over  the  EVT  threshold  (depending  on  the  SNR).  These  results  were  verified  with  a  simulated  data 
sequence  and  real  fade-rate  data  obtained  over  a  16.2  km  link  (32.4  km  round  trip  off  of  retrorefiectors)  at  the  NRL  FSO 
facility.  These  preliminary  findings  indicate  that  adaptive  thresholding  may  be  implemented  to  improve  BER  in  a  free 
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space  optical  link  with  significantly  less  intensive  computational  loads  and  commensurately  lower  costs.  In  future  work, 
we  will  conduct  experiments  to  verify  these  results. 
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